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Abstract. Recent transport experiments in the original cuprate high temperature 
superconductor, La2-xBa2,Cu04, have revealed a remarkable sequence of transitions 
and crossovers which give rise to a form of dynamical dimensional reduction, in which 
a bulk crystal becomes essentially superconducting in two directions while it remains 
poorly metallic in the third. We identify these phenomena as arising from a distinct new 
superconducting state, the "striped superconductor," in which the superconducting 
order is spatially modulated, so that its volume average value is zero. Here, in addition 
to outlining the salient experimental findings, we sketch the order parameter theory 
of the state, stressing some of the ways in which a striped superconductor differs 
fundamentally from an ordinary (uniform) superconductor, especially concerning its 
response to quenched randomness. We also present the results of DMRG calculations 
on a model of interacting electrons in which sign oscillations of the superconducting 
order are established. Finally, we speculate concerning the relevance of this state to 
experiments in other cuprates, including recent optical studies of La2-a;Sr2;Cu04 in a 
magnetic field, neutron scattering experiments in underdoped YBa2Cu306+x, and a 
host of anomalies seen in STM and ARPES studies of Bi2Sr2CaCu208+5- 
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1. Introduction 

In this paper we carefully characterize in terms of its broken symmetries, a novel 
superconducting state of matter, the "pair-density- wave" (PDW), with special focus on 
the "striped superconductor" , a unidirectional PDW. We present a concrete microscopic 
model of interacting electrons which we show, using density matrix renormalization 
group (DMRG) methods, has a striped superconducting ground state. There is an 
intimate relation between PDW and charge density wave (CDW) order, as a consequence 
of which the striped superconductor exhibits an extreme sensitivity to quenched disorder 
which inevitably leads to glassy behavior. This is qualitatively different from the familiar 
effects of disorder in uniform superconductors. 

On the experimental side, we first draw attention to a set of recently discovered 
transport anomahes in the high temperature superconductor, La2-a;Baa;Cu04, which 
are particularly prominent for x = 1/8. We will be particularly interested in the 
spectacular dynamical layer decoupling effects recently observed in this system [1] which 
indicate that the effective inter-layer Josephson coupling becomes vanishingly small with 
decreasing temperature. These experiments suggest that a special symmetry of the state 
is required to explain this previously unsuspected behavior. While no comprehensive 
theory of these observations currently exists, even at the phenomenological level, we 
show how the salient features of these observations can be straightforwardly rationalized 
under the assumption that La2_a;Ba2:Cu04 is a striped superconductor. We outline 
some further experiments that could critically test this assumption. Finally, we 
speculate about the possible role of striped superconducting order as the source of 
a number of salient experimental anomalies in a much broader spectrum of high 
temperature superconductors, including recent experiments on magnetic field induced 
layer decoupling in La2_2:Sr^Cu04 [2], the notable evidence of a local gap with 
many characteristics of a superconducting gap in STM and ARPES experiments in 
Bi2Sr2CaCu208+5 [3llll|5], and, most speculatively of all, experiments indicative of time 
reversal symmetry breaking in the pseudo-gap regime of YBa2Cu306+a;. 

The striped superconductor is a novel state of strongly correlated electronic matter 
in which the superconducting, charge and spin order parameters are closely intertwined 
with each other, rather than merely coexisting or competing. As we show here (and 
discussed in [6l [7j) the striped superconductor arises from the competing tendencies 
existing in a strongly correlated system, resulting in an inhomogeneous state in which 
all three forms of order are simultaneously present. The striped superconductor is 
thus a new type of electronic liquid crystal state [Sj.f In particular, as we shall 



f Electronic liquid crystals [8] are quantum states of matter that spontaneously break some of the 
translation and/or rotational symmetries of an electronic system. In practice, these symmetries are 
typically not the continuous symmetries of free space, but rather the various discrete symmetries of the 
host crystal. Although an electronic smectic (stripe ordered state) has the same order parameter as a 
CDW or an SDW, it is a more general state that does not necessarily derive from a nesting vector of 
an underlying Fermi surface. The liquid crystal picture offers a broader perspective on the individual 
phases and on their phase transitions [3] ■ In particular, the way in which the PDW phase intertwines 
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see, the symmetry breaking pattern of the striped superconductor naturally explains 
the spectacular layer decoupling effects observed in experiments in La2_xBaj;Cu04. 
In contrast, in any state with uniform superconducting order, dynamical inter-layer 
decoupling could only arise if a somewhat unnatural fine tuning of the inter-layer 
couplings led to a sliding phase fiO\ . 

The striped superconductor has an order parameter describing a paired state with 
non- vanishing wave vector, Q, the ordering wave vector of the unidirectional PDW. 
As such, this state is closely related by symmetry to the Fulde-Ferrell [H] (FF) and 
Larkin-Ovchinnikov [12] (LO) states. The order parameter structure of the PDW state, 
involving several order parameters coupled to each other, also evokes the 50(5) approach 
of a unified description of antiferromagnetism and a uniform rf-wave superconductor 
[T3l [Hj. The relation of the present discussion to these other systems and to earlier 
theoretical works on the same and closely related subjects is deferred to Sec. El The 
physics of stripe phases in the cuprate superconductors has been reviewed in Ref.[T5] 
and more recently in Ref. [16]. 

It is important to stress that the macroscopically superconducting phase of the 
cuprates reflects the existence of a spatially uniform Q = component of the order 
parameter, whether or not there is substantial finite range superconducting order at 
non-zero Q. One might therefore reasonably ask whether striped superconductivity, 
even if interesting in its own right, is anything but an exotic oddity, with little or no 
relation to the essential physics of high temperature superconductivity. The answer to 
this question is at present unclear, and will not be addressed to any great extent in the 
present paper. However, we wish to briefly speculate on a way in which the striped 
superconductor could play an essential role in the broader features of this problem. In 
BCS theory, the superconducting state emerges from a Fermi liquid in which the strong 
electronic interactions are already accounted for in the self-energy of the quasiparticles. 
The cuprates are different, in that the superconductivity, especially in underdoped 
materials, emerges from a pseudogap phase for which there is no commonly accepted 
model. As we will show, striped superconductivity has features in common with the 
pseudogap phase, such as a gapless nodal arc and antinodal gaps. We speculate that the 
pseudogap phase might be associated with fluctuating striped superconductivity, a state 
that we do not yet know how to treat. Nevertheless, analysis of the ordered PDW state 
and comparison to observations of stripe-ordered cuprates is a starting point. Indeed, 
comparison (see Subsection 19.11) of recent photoemission results on LBCO x=l/8 with 
transport and optical properties suggests that a "uniform" d-wave state {i.e. one with 
a non-zero uniform component of the order parameter) develops on top of a striped 
superconductor, resulting in a fully superconducting Meissner state, albeit one with 
substantial coexisting short-range correlated stripe order. 

The rest of this paper is organized as follows: In Section [2] we give an order 
parameter description of the pair-density wave state. In Sees. [3] and HI we summarize 

charge, spin, and superconducting orders is unnatural in terms of a Fermi surface instability, but not 
so from the liquid crystalline perspective. 
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the experimental evidence for this state, with Sec. [3] focussing on the strongest case, 
La2_2:Baa;Cu04, and Sec. [Hon other cuprates. In Section [5] we discuss the microscopic 
mechanisms for the formation of a PDW state. In Subsection 15.11 we implement these 
microscopic considerations by constructing a specific model that exhibits a PDW phase. 
The central conceptual ingredient is a microscopic mechanism leading to the formation 
of TT junctions in an unidirectional PDW state, which is given in Subsection 15.1.11 
using perturbative arguments and then checked numerically using the density matrix 
renormalization group (DMRG) (in Subsection 15.1.21) . A solvable microscopic model is 
discussed in Subsection 15. 2[ The quasi-particle spectrum of the PDW state is discussed 
in Subsection 15.31 Next, the Landau-Ginzburg theory of the PDW phase is discussed 
in Section [61 In Section [7] we show that the PDW state, in three dimensional layered 
structures (orthorhombic and LTT) as well as at grain (twin) boundaries, leads to 
time-reversal symmetry breaking effects. In Section [8] we discuss the connections that 
exist between the PDW state and other states discussed in the literature, particularly 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states. Section [9] is devoted to our 
conclusions. 

This paper is partly a review of our recent work on the theory of the PDW state[6l[7] 
and of other related work, with an updated discussion of the current experimental 
status. However in this paper we have also included many new results, particularly the 
DMRG analysis of PDW states in strongly correlated systems of Section 15.1.21 and the 
connection bewteen the PDW state and non-coUinear order and time-reversal symmetry 
breaking of Section [71 

2. The Order Parameter of a Striped Superconductor 

The order parameter whose non-zero expectation value defines a superconducting state 
is 



where ipl{r) is the fermionic field operator which creates an electron with spin 
polarization a at position r. Further distinctions between different superconducting 
states can be drawn on the basis of the spatial and spin symmetries of 0. In crystalline 
solids, all familiar superconducting states respect the translational symmetry of the 
solid, 0(r + R, r' + R) = 0(r, r'), where R is any Bravais lattice vector. Consequently, 
the symmetries of the state can be classified by the irreducible representations of the 
point group - colloquially as s-wave, d-wave, p-wave, etc. In the absence of spin-orbit 
coupling, superconducting states can be classified, as well, by their transformation under 
spin rotations as singlet or triplet. Finally, the superconducting state can either preserve 
or break time reversal symmetry (as in p^ + ipy)- 

In the presence of quenched disorder, the underlying Hamiltonian does not have any 
particular spatial symmetries, so the classification of distinct superconducting states by 
their symmetries (other than time reversal), at first seems difficult. However, there are 




(1) 



Striped superconductors 



5 



several ways that this can be accomphshed [T7j, of which the most obvious is to consider 
the symmetries of the configuration averaged order parameter 

0,,,,(r,r') = (4(r)V^Mr')), (2) 

where (...) signifies the thermal average, and (. . .) signifies an average over realizations 
of the disorder configuration. It is clear, for example, that under most circumstances, 
a macroscopic "phase sensitive" measurement of the symmetry of the order parameter 
will give [T7] a result consistent with a classification based on the symmetry of the 
configuration averaged order parameter. 

The striped superconductor is an example of a state, which has more generally 
called [TSl [T9l [20] a pair density wave (PDW), in which the translational symmetry of 
the crystal is spontaneously broken as well, so that (j){r + R, r' + R) exhibits non-trivial 
dependence on R. However, this by itself, is insufficient to identify a new state of matter. 
In a system with coexisting charge- density- wave (CDW) and superconducting order, the 
CDW itself introduces a new periodicity into the problem, which must generically be 
reflected in a spatial modulation of 0, as well.f As discussed in Sec. El an analysis of 
the implications of a generic theory of coupled order parameters implies [7] that in a 
state of coexisting order, a (possibly small) modulation of the superconducting order 
with the same spatial period as that of the CDW will be induced. None-the-less, in 
such a state, there still exists a "dominant" uniform component to the superconducting 
order parameter, which we define as the spatial average of the SC order parameter: 

0S;i,(r, r') ^ iV-i ^(V^t (r + K)^l,{r' + R)), (3) 

R 

where N is the number of unit cells in the system. 

Instead, the pure PDW in a crystal is a state in which is non-zero, but all uniform 
components vanish, 0^°^, (r,r') = for any r and r'. Just as a CDW is often defined 
in terms of a fundamental harmonic, so a PDW state is characterized by the smallest 
value of the crystal momentum, Q, for which 

0i5(r, r') = J2 e^P[^Q ■ ^Mi^ + R)V^l'(r' + R)), (4) 

R 

has a non- vanishing expectation value. | 

f The problem of coexisting stripe and superconducting order in strongly correlated systems has been 
the focus of numerous studies in the literature. Sachdev, Vojta, and coworkers have investigated this 
problem in detail in the context of generalized 2D t — J models in the large N approximation [2T 1 \22 \ I23j . 
This problem has also been discussed in one-dimensional systems [24]. 

I As with a uniform superconducting state, distinct PDW states with the same pattern of translation 
symmetry breaking can also be distinguished by different patterns of point group symmetry breaking. 
However, since the ordering vector (or vectors) already break the point group down to a smaller 
subgroup, which is then all that is left of the original symmetry for this purpose. For instance, in a 
tetragonal crystal, a striped superconductor with Q along the x direction, can be classified as having 
s-wave or da;j,-wave symmetry, based on whether or not the the order parameter changes sign under 
reflection through a symmetry plane parallel to the x axis, but any distinction one would like to draw 
between a striped version of an s-wave and a d2;2_j,2-wave superconductor are in precise, not based on 
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Note that the theory of coupled order parameters [7j imphes that the existence of 
PDW order with ordering vector Q generically imphes the existence of CDW order with 
ordering vector 2Q, but so long as 0° = 0, no CDW ordering with wave vector Q is 
expected. A "striped superconductor" refers to the special case in which the independent 
ordering vectors are all parallel to each other ("unidirectional PDW"). 

One of the prime new characteristics of a striped superconductor which is different 
from a uniform superconductor is its complex sensitivity to quenched disorder. As we 
shall see, for much the same reasons that disorder destroys long-range CDW order, under 
most relevant circumstances, even weak disorder causes the configuration averaged PDW 
order parameter to vanish: 

0.,.'(r,r') = O. (5) 

However, as in the case of an XY spin-glass, this is not the whole story: It is possible 
to define an analogue of the Edwards- Anderson order parameter, 

g.,.,(r,r') = |(4(r)^l,(r'))P, (6) 

which vanishes in the normal high temperature phase, but which can be non-zero in a 
low temperature superconducting glass phase, where one exists. Moreover, in such a 
phase, as we will see, we generically expect time-reversal symmetry to be spontaneously 
broken. In analogy with the XY spin-glass, we expect that in two dimensions, the 
superconducting glass phase is stable only at T = and for weak enough disorder, 
although in three dimensions it can exist below a non-zero superconducting glass 
transition temperature. § 

There is one more extension that is useful — we define a charge 4e superconducting 
order parameter: 

0(^)(1,2,3,4) ^ {^Ur,)<i^2)€A^3)^U^4)) (7) 

where we have introduced a compact notation in which 1 = (cri,ri), etc. Naturally, in 
any state with charge 2e superconducting order, (j)„y{r,r') ^ 0, some components of 
the charge 4e order parameter will also be non-zero. This can be seen from the theory 
of coupled order parameters presented in Sec. [61 At mean-field level, it can be seen by 
applying Wick's theorem to the expression in Eq. [7] to express 0'^^'' as a sum of pairwise 
products of 0's: 2, 3, 4) ~ 0(12)0(34) + 0(14)0(23) - 0(13)0(24). 

There are two reasons to consider this order parameter. In the first place, it is 
clear from the above that even in the PDW state, although the uniform component 
of vanishes, the uniform component of 0^ ~ 0^0^*-^ ^ 0. More importantly, 
0(4) ordering can be more robust than the PDW ordering. Specifically, under some 

broken symmetries but on quantitative differences in local pairing correlations, and so do not define 
distinct phases of matter. However, a checkerboard PDW, with symmetry related ordering vectors Q 
and Q' along the x and y axes, respectively, can be classified as s-wave or (ia.2_j^2-wave, depending on 
how it transforms under rotation by 7r/2. 

§ In [7J, the possibility is discussed that in three dimensions there might also exist a superconducting 
version of a Bragg glass phase, in which exhibits quasi-long-range order. We have not further studied 
this potentially interesting state. 
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circumstances [71 [2S], it is possible for thermal or quantum fluctuations to destroy the 
PDW order by restoring translational symmetry without restoring large-scale gauge 
symmetry; in this case, appropriate components of (j)^^^ remain non-zero, even though 
vanishes identically. This is the only potentially realistic route we know of to charge 4e 
superconductivity. II 

In the absence of spin-orbit coupling, distinct phases with translationally invariant 
charge 4e ordering can be classified according to the total spin of the order parameter, 
which in this case can be spin 2, 1, or 0. Manifestly, any charge 4e superconducting 
state which results from the partial melting of a singlet PDW will itself have spin 0. 
As with paired superconductors, the charge 4e order parameter can also be classified 
according to its transformation properties under action of the point-group of the host 
crystal. For instance, in a crystal with a C4 symmetry, taking the points to lie on the 
vertices of a square, the transformation properties of 0*^^^ under rotation by 7r/2 can be 
used to classify distinct spin-0 states as being d-wave or s-wave. 

The definitions given here in terms of possible behaviors of the order parameter are 
natural from a taxonomic viewpoint. In particular, the striped superconductor seems at 
first to be a rather straightforward generalization of familiar uniform superconducting 
states. However, both at the microscopic level of the "mechanism" of formation of such 
a state, and at the phenomenological level of macroscopically observable implications of 
the state, the problem is full of subtleties and surprises, as discussed below. 

3. Striped superconductivity in La2 2;Ba^Cu04 and the 214 family 

We now summarize some of the observations that lead to the conclusion that 
La2-2:Ba2:Cu04 with x = 1/8 is currently the most promising candidate experimental 
system as a realization of a striped superconductor [1], [6] . 

Firstly, the existence of "stripe order" is unambiguous. It is well known (from 
neutron [271 EHl EH] and X-ray [301 ED EH] scattering studies) that unidirectional 
CDW (charge-stripe) and SDW (spin-stripe) orders exist in La2_a;Baa;Cu04. Such 
spin and charge stripe orders were originally studied in Lai,48Ndo.4Sro.i2Cu04 [32l 
l33l and they have now been confirmed in Lai,8-a;Euo.2Sra,Cu04 [35], [361 137] . 
Furthermore, substantial spin stripe order has been observed in La2CuO4.11 [38], Zn- 
doped La2-a;Sra;Cu04 [39J], and in the spin-glass regime of La2-a;Sra;Cu04 (0.02 < 
X < 0.055) [10]. While the spin-stripe order in underdoped but superconducting 
La2-2:Sr2;Cu04 (0.055 < a; < 0.14) is weak in the absence of an applied magnetic 
field, it has been observed [HI HJl |13] that readily accessible magnetic fields (which 
partially suppress the superconducting order) produce well-developed and reproducible 
spin-stripe order. For La2-xBaa.Cu04 with x = 1/8, the charge ordering temperature is 

II It is possible to cook up models in which charge 4e superconductivity arises in systems in which 
electrons can form 4 particle bound-states, but do not form 2 particle or many particle bound states 
(phase separation) - see, for example, However, this involves unrealistically strong attractive 

interactions and an unpleasant amount of fine tuning of parameters. 
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53 K and the spin ordering temperature is 40 K |28j . 

In addition to spin and charge ordering, La2_a,.Ba2:Cu04 with x = 1/8 exhibits 
transport and thermodynamic behavior that is both striking and complex. We will not 
rehash all of the details here (see ^28^); however, there are two qualitative features of the 
data on which we would like to focus: 1) With the onset of spin-stripe order at 40 K,f 
there is a large (in magnitude) and strongly temperature dependent enhancement of 
the anisotropy of the resistivity and other properties, such that below 40 K the in- 
plane charge dynamics resemble those of a superconductor, while in the c-direction 
the system remains poorly metallic. The most extreme illustration of this occurs 
in the temperature range 10 K < T < 16 K, in which the in-plane resistivity is 
immeasurably small, while the c-axis resistivity is in the 1-10 mfi-cm range, so that 
the resistivity anisotropy ratio is consistent with infinity. 2) Despite the fact that 
signatures of superconductivity onset at temperatures in excess of 40 K, and that angle 
resolved photoemission (ARPES) has inferred a "gap" [HlllS] of order 20 meV, the fully 
superconducting state {i.e., the Meissner effect and zero resistance in all directions) only 
occurs below a critical temperature of 4 K. It is very difficult to imagine a scenario in 
which a strong conventional superconducting order develops locally on such high scales, 
but fully orders only at such low temperatures in a system that is three dimensional, 
non-granular in structure, and not subjected to an external magnetic field. 

Evidence that similar, although somewhat less extreme transport and thermody- 
namic anomalies accompany stripe ordering can be recognized, in retrospect, in other 
materials in the 214 family. For example, in La2-a:Sra;Cu04 with x = 0.08 and 0.10, 
the anisotropy of the resistivity (c-axis vs. in-plane) rapidly grows towards lO'^ as the 
superconducting Tc is approached from above [46j. In the case of Lai.6-xNdo.4Sra;Cu04, 
evidence for dynamical layer decoupling is provided by measurements of the anisotropic 
onset of the Meissner effect [37]. In contrast, the resistivity ratio in this material [H] only 
reaches 10^; this may be limited by enhanced in-plane resistivity due to disorder |49j . 
Moreover, an unexpectedly strong layer decoupling in the charge dynamics produced by 
the application of a transverse magnetic field in La2-a;Sra;Cu04 has been observed [2], 
but only in the underdoped range of x where the magnetic field also induces spin stripe 
order [4ll|43]. 

We shall see that the anomalous sensitivity of a striped superconductor to quenched 
disorder can account for the existence of a broad range of temperatures between 
the onset of strongly developed superconducting correlations on intermediate scales 
and the actual macroscopic transition temperature to a state of long-range coherence. 
Moreover, given the crystal structure of the Low Temperature Tetragonal (LTT) phase of 
La2_2:Ba2:Cu04, there is a special symmetry of the striped superconducting state which 
produces interlayer decoupling. Specifically, because the stripes in alternate planes 
are oriented perpendicular to one another (as shown in FigJS^), there is no first order 
Josephson coupling between neighboring planes. Indeed, analogous features of the spin- 

t Static charge stripe order onsets at 54 K, at the LTT/LTO transition. [30] 
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stripe order, which have gone largely unnoticed in the past, are accounted for [B] by 
the same geometric features of the striped state. When spin-stripe order occurs, the 
in-plane correlation length can be very long, in the range of 100-600 A [381 EH], but 
the interplanar correlation length is never more than a few A [38l [50] , a degree of 
anisotropy that cannot be reasonably explained simply on the basis of the anisotropy 
in the magnitude of the exchange couplings [51] . Furthermore, despite the presence of 
long correlation lengths, true long-range spin-stripe order has never been reported. It 
will be made clear that the superconducting stripe order and the spin stripe order share 
the same periodicity and the same geometry, so both the interlayer decoupling and the 
suppression by quenched disorder of the transition to a long-range ordered state can be 
understood as arising from the same considerations applied to both the unidirectional 
SDW and PDW orders. 

4. Experiments in other cuprates 

Data which clearly reveal the existence of spin-stripe order, or that provide compelling 
evidence of of PDW order in other families of cuprates is less extensive. However, 
considerable evidence of a tendency to spin stripe order in YBa2Cu306+x has started to 
accumulate |15j. and there are some persistent puzzles concerning the interpretation of 
various experiments in a number of cuprates that, we would like to speculate, may reflect 
the presence of PDW order. In this section, we will mention some of these puzzles, and 
will return to discuss why they may be indicative of PDW order in Sec. [71 

YBa2Cu306+a: is often regarded as the most ideal cuprate, having minimal 
structural and chemical disorder, and less tendency to stripe or any other type of charge 
ordering than the 214 cuprates. (Sometimes the cuprates with the highest transition 
temperatures, such as IIgBa2Cu3 044.5, are viewed as being similarly pristine.) However, 
in its underdoped regime it is well known that YBa2Cu306+a; exhibits temperature- 
dependent in-plane anisotropic transport [52J as well as fluctuating spin stripe order 
[531 1151 Recent neutron scattering experiments have provided strong evidence 

that underdoped YBa2Cu306+a; (with x ~ 0.45) has nematic order below a critical 
temperature ~ 150 K ^55j. Even more recent neutron scattering experiments 
by Hinkov et al. [56] on the same sample find that a modest c-axis magnetic field 
stabihzes an incommensurate static spin ordered state, detectable as a pair of peaks 
in the elastic scattering displaced by a distance in the crystallographic a direction 
from the Neel ordering vector. Given the newfound evidence of spin-stripe related 
structures in YBa2Cu306+x, it is plausible that here, too, striped superconductivity 
may occur. However, the differences in the 3D crystal structure, and especially the 
weak orthorhombicity, would make the macroscopic properties of a PDW distinctly 
different in YBa2Cu306+2; than in the 214 cuprates. 

A remarkable recent discovery is that underdoped YBa2Cu306+z appears to exhibit 
signatures of spontaneous time-reversal symmetry breaking (at zero magnetic field) 
below a critical temperature comparable to that for the nematic ordering. [571 EH] 
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(HgBa2Cu304+5 [SS] exhibits similar signatures.) Various theoretical scenarios for the 
existence of time-reversal symmetry breaking predated these experiments, and so in 
some sense predicted them [601 [6T] . However, given that both nematic order and 
time-reversal symmetry breaking are seemingly present simultaneously in the same 
samples with comparable critical temperatures, it is reasonable to hope that both 
phenomena have an underlying common explanation. If we think of the superconducting 
order parameter as an XY pseudo-spin, then the PDW order is a form of coUinear 
antiferromagnetism, and time-reversal symmetry breaking corresponds to non-coUinear 
order of the pseudospins. As we will show in Sec. [71 weak time reversal symmetry 
breaking can occur in a PDW state due to various patterns of geometric frustration 
in three dimensions or as a consequence of the existence of certain types of defects, 
such as twin boundaries. (See, also, [7]. ) There is a large body of STM and 
ARPES data, especially on Bi2Sr2CaCu208+6 and Bi2Sr2Cu06+5, which has revealed 
a surprisingly rich and difficult to interpret set of spectral features associated with the 
d-wave superconducting gap and a d-wave pseudo-gap whose origin is controversial. 
Indeed, there is a clear "nodal-anti nodal dichotomy" [62l [63] in the behavior of the 
measured single-particle spectral functions. Some aspects of the data are suggestive that 
there is a single superconducting origin of all gap features, with anisotropic effects of 
superconducting fluctuations leading to the observed dichotomy. Other aspects suggest 
that there are at least two distinct origins of the near-nodal and the antinodal gaps. It is 
possible that PDW ordering tendencies can synthesize both aspects of the interpretation. 
In the presence of both uniform and PDW superconducting order, there are two distinct 
order parameters, both of which open gaps on portions of the Fermi surface, but they 
are both superconducting, and so they can smoothly evolve into one another. (Note 
that an early study [64J of modulated structures seen in STM [65|, [3l [66] concluded that 
they could be understood in terms of just such a two-superconducting-gap state.) 

More generally, one of the most remarkable features of the pseudo-gap phenomena is 
the existence of what appears to be superconducting fluctuations, detectable [671 [Ml [69] 
for instance in the Nernst and magnetization signal, over a surprisingly broad range of 
temperatures and doping concentrations. At a broad-brush level [70] , these phenomena 
are a consequence of a phase stiffness scale that is small compared to the pairing scale. 
However, it is generally difficult to understand the existence of such a broad fiuctuational 
regime on the basis of any sensible microscopic considerations. The glassy nature of 
the ordering phenomena in a PDW may hold the key to this central paradox of HTC 
phenomenology, as it gives rise to an intrinsically broad regime in which superconducting 
correlations extend over large, but not infinite distances. 

5. Microscopic considerations 

From a microscopic viewpoint, the notion that a PDW phase could be stable at first 
sounds absurd. Intuitively, the superconducting state can be thought of as the condensed 
state of charge 2e bosons. However, in the absence of magnetic fields, the ground-state 
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of a bosonic fluid is always node-less, independent of the strength of the interactions, 
and therefore cannot support a state in which the superconducting order parameter 
changes sign. Thus, for a PDW state to arise, microscopic physics at scales less than or 
of order the pair-size, must be essential. This physics reflects an essential difference 
between superfluids of paired fermions and preformed bosons |71] . 

Our goal in this section is to shed some light on the mechanism by which strongly 
interacting electrons can form a superconducting ground-state with alternating signs 
of the order parameter. We will consider the case of a unidirectional (striped) 
superconductor, but the same considerations apply to more general forms of PDW 
order. We will not discuss the origin of the pairing which leads to superconductivity. 
Likewise, we will not focus on the mechanism of translation symmetry breaking by the 
density wave, as that is similar to the physics of CDW and SDW formation. Our focus 
is on the sign alternation of 0. Thus, in much of this discussion, we will adopt a model 
in which we have alternating stripes of superconductor and correlated insulator. The 
system looks like an array of extended superconductor-insulator-superconductor (SIS) 
junctions, and we will primarily be concerned with computing the Josephson coupling 
across the insulating barriers. If the effective Josephson coupling is positive, then a 
uniform phase (normal) superconducting state is favored, but if the coupling is negative 
(favoring a vr junction), then a striped superconducting phase is found. 

So long as time reversal symmetry is neither spontaneously nor explicitly broken, 
the Josephson coupling, J between two superconductors must be real. If it is positive, 
as is the usual case, the energy is minimized by the state in which the phase difference 
across the junction is 0; if it is negative, a phase difference of vr is preferred, leading to 
a "vr junction." vr junctions have been shown, both theoretically and experimentally, to 
occur for two distinct reasons: they can be a consequence of strong correlation effects 
in the junction region between two superconductors [TTl [72| [73] or due to the internal 
structure {e.g. d-wave symmetry) of the superconductors, themselves [TillTS]. 

In Ref. [7], we have provided examples of vr junctions which build on the first set 
of ideas. Unlike the previously studied cases, these tt junctions were extended (i.e., J is 
proportional to the cross sectional "area" of the junction). However, since the problem 
was solved analytically (treating the tunneling between the superconducting and the 
insulating regions by perturbation theory), we were limited to somewhat artificial 
models. For example, tunneling between the sites in the insulating regime was neglected. 
In this Section, we first summarize the perturbative results of [7], and then present 
numerical (DMRG) results for an extended SIS junction. Under some circumstances, 
J > 0, but we also find a considerable region of parameter space where where J < 0. 
Finally, we discuss how this result can be generalized to an infinite array of junctions, 
forming a 2D unidirectional PDW. 
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5.1. A Solved model 

Let us consider the following explicit model for a single SIS junction. The three 
decoupled subsystems are described by the Hamiltonian 

Ho = Hl + Hb + Hr, (8) 

The right {R) and left (L) superconducting regions and the barrier (B) region are one 
dimensional Hubbard models, 

i(7 i 

i+i cr ^ creation operator of an electron on chain a = L, i? or i? at site i with spin a, 
and we have introduced the notation na,i,o- = '^a,i,a'^a,i,(j ^'^d na,i = Ylia^a,i,a- The left 
and right superconducting chains are characterized by a negative Ur = Ul = — \ Ul^r\, 
while the insulating barrier has a positive f/^ > 0. The chemical potentials of the left 
and right superconductors are the same, iir = fii, but different from /i^, which is tuned 
such that the barrier chain is half filled (and therefore insulating). 

The three subsystems are coupled together by a single-particle hopping term, 

H' = -t' Et^l.'^^B.,. + C^^,aCB,^,. + h-C]- (10) 
i,(7 

The left and right chains are characterized by a spin gap and by dominant 
superconducting fluctuations, as a result of their negative t/'s. The inter-chain hopping 
term H' induces a finite Josephson coupling between the local superconducting order 
parameters of the two chains, via virtual hopping of a Cooper pair through the barrier 
chain. 



5.1.1. Perturbative analysis For completeness, let us briefly review the perturbative 
treatment of the inter-chain hopping term (fTOj) given in [7]. The leading (fourth order) 
contribution to the Josephson coupling is given by 

it')' 



J 



f3 



where 1 = {Ti,ii) etc., 
(in the limit (3 oo) and 



dl d2 d3 dA Fl{1, 2)F^(4, 3)1(1, 2; 3, 4) 



dri 



Fa(l,2) 



r(l,2;4,3) 



Tr 
Tr 



cl,n,t(n)c^a(^2) 

4i,T(^i)4,i(^2)c.3,|(r3)Q4,t(r4) 



(11) 



(12) 



(13) 



where we have made the identification c|,^ = c^j^.. Our purpose is to determine the 
conditions under which J < 0. For the sake of simplicity, let us consider the case in 
which the gap to remove a particle from the barrier, A/^, satisfies <^ A/i <C Ap, 
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where is the spin gap on the superconducting chains, and Ap is the gap to insert 
a particle in the barrier. These conditions can be met by tuning appropriately the 
chemical potentials on the three chains and setting Ub to be sufficiently large. 

In [7] it is shown that, quite generally, J can be written as a sum of two terms 

J = Jl + J2 (14) 

— * — # 

where, in terms of the spin-spin correlation function, (5(1) ■ S{2)) of the barrier chain. 



it 



f\4 



Ji = ^" ' n I dl d2\FL{l,2)? (15) 
4/3(A,)2y ' ' ^' ^ ' 

J2=--|^/ dl d2\F,{l,2)\^S{l)-S{2)) 

Explicitly, Ji > 0, while for generic circumstances one finds that J2 < 0. The overall 
sign of J is therefore non-universal, and determined by which term is bigger. We can, 
however, identify the conditions under which J2 dominates. Upon a Fourier transform, 
|Fi(l, 2)1^ is peaked around two values of the momentum g, at g = and 2fcj?, in which 
2/ci? = vrn where n is the number of electrons per site in the left and right chains. Since, 
upon Fourier transforming, (5(1) ■ 5(2)) is peaked at momenta g = and vr, (as can be 
seen, e.g., from a bosonized treatment of the half ffiled chain) we expect that J2 in Eq. 
f fTSj) is maximized when n = 1, i.e. when the superconducting chains are half filled. The 
requirement of proximity to half ffiling becomes less and less stringent when |?7ij,L| is 
increased, since then the gap A^ in the superconducting chains increases and the peaks 
in \FL[q)\^ become more and more broad. These qualitative expectations are conffimed 
by numerical DMRG simulations, presented in the next subsection. 



5.1.2. Numerical results We have performed DMRG simulations of the model H = 
Hq + H' in Eq. flHfTOj) . with the following parameters: t = t' = l,/is = 6, [7^ = 10, and 
variable Ul = Ur = — \ Ul,r\ and fiL = fiR = fiL,R- Most of the calculations where done 
with systems of size 3 x 24. In a small number of parameter sets, we have verified that 
the results do not change when we increase the system size to 3 x 36. Up to m = 1600 
states where kept in these calculations. The results (both ground state energies and 
local measurements) where extrapolated linearly in the truncation error [76], which is 
in the range 10~^ — 10~^. 

In order to measure the sign of the Josephson coupling from the calculations, we 
have applied pairing potentials on the left and right chains, adding the following term 
toEq. dHD: 

Hpair = - Yl ^c.^i,AM + '^-^- (16) 

i,a=L,R 

In the presence of this term, the number of particles in the calculation is conserved only 
modulo 2. The average particle number is fixed by the overall chemical potential. Two 
methods where employed to determine the sign of J. (a) Pairing potentials of either 
the same sign, Ar = A^, and of opposite signs, Ar = —Ai, where applied to the two 
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chains. The ground state energies in the two cases are and E_, respectively. Then 
J = — [77]. (b) A pairing potential was applied to the left chain only, > 0, 
while Ar = 0. The induced pair field 



on the right chain was measured. Its sign indicates the sign of J. This is the method 
we used in most calculations. Method (a) was applied to a small number of points in 
parameter space, and found to produce identical results to those of method (b) for the 
sign of J. 

Fig. [T] shows the local expectation values of the particle number, spin and pair 
field operators along the three chains for \Ul,r\ = 2.5 and various values of ^lr- The 
density of electrons on the left and right chains increases as fi^R increases, while the 
density on the middle chain is kept close to one particle per site. A positive pair 
potential of strength = 0.1 was applied on the left chain, inducing a positive pair 
field 0L = (c|^j|C^j|) > 0, while Ar = 0. A negative induced pair field (j)R on the right 
chain indicates that the effective Josephson coupling J between the left and right chains 
is negative. Note that J is negative for the two upper rows (in which (^l,/?) = 0.9, 
0.83 respectively), while for the two lower rows (where (^l,/?) = 0.75, 0.66) it becomes 
positive. This is in agreement with our expectation, based on the perturbative analysis 
of the previous subsection, that when the superconducting chains are close to half filling 
{{n) = 1), the negative J2 term dominates and the overall Josephson coupling is more 
likely to become negative. 

The middle column in Fig. [1] shows the expectation value of the z component of 
the spin along the three chains. In order to visualize the spin correlations, a Zeeman 
field of strength h = 0.5 was applied to the i = 1 site of the middle chain. The results 
clearly indicate that the two outer chains have a spin gap (and therefore have a very small 
induced moment), while in the half filled middle chain there are strong antiferromagnetic 
correlations. Interestingly, as the chemical potential on the outer chains is decreased, the 
spin correlations along the middle chain become incommensurate. This seems to occur 
at the same point where the Josephson coupling changes sign (between the second and 
third row in Fig. [1]). This phenomenon was observed for other values of Ul^r, as well. 
The incommensurate correlations can be explained by the further-neighbor Ruderman- 
Kittel-Kasuya-Yosida (RKKY)-like interaction which are induced in the middle chain 
by the proximity of the outer chains. Upon decreasing the inter-chain hopping t' to 0.7, 
the spin correlations in the middle chain become commensurate over the entire range 
of fiL,R (and the region of negative J increases). Why J > seems to be favored by 
incommensurate correlations in the middle chain is not clear at present. 

Fig. [2] shows the phase diagram of the three chain model as a function of the density 
{fiLfi) and the attractive interaction \Ul,r\ on the outer chains. In agreement with the 
perturbative considerations, proximity to {til r) = 1 and large \Ul^r\ (compared to the 
bandwidth 4t) both favor a negative Josephson coupling between the outer chains. 




(17) 
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Figure 1. (Color online.) The left, middle and right columns show the average 
density {rii), z component of the spin (Sf) and pair field (ci^|Ci_|), respectively, as 
a function of position i along the chains, calculated by DMRG for 3 x 24 systems. 
Circles, diamonds and dots refer to the left, middle and right chains, respectively. The 
attractive interaction on the superconducting (left and right) chains is \Uh.B\ — 2.5 in 
all calculations. A pairing term [Eq. (fTB|) ] was applied with = 0.1 and Aj^ — 0. 
The other model parameters are given by: i = i' = 1, /xs = 6, t/s = 10. Each row 
corresponds to a single calculation with a specific value of the chemical potential /xl,^ 
(and hence a particular particle density) on the superconducting chains. 



5.2. Extension to an infinite array of coupled chains 

The model presented in the previous subsection includes only a single extended tt 
junction. However, it is straightforward to extend this model to an infinite number 
of coupled chains with alternating U. So long as the Josephson coupling across a single 
junction is small, we expect that the extension to an infinite number of chains will 
not change it by much. Therefore, in the appropriate parameter regime in Fig. [2l the 
superconducting order parameter changes sign from one superconducting chain to the 
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Figure 2. (Color online.) Phase diagram of the three chain model [Eq. (|8ll0p ] from 
DMRG, as a function of \Ul.b\j the attraction on the left and right chains, and {tll^h), 
the number of electrons per site on the left and right chains. The following parameters 
were used: t = t' — 1, — 6, Ub = 10. On the middle chain, (n^) w 1 in all cases. 
The symbols show the points that were simulated. 

next, forming a striped superconductor (or unidirectional PDW). 

In order to demonstrate that there are no surprises in going from three chains to 
two dimensions, we have performed a simulation for a 5 x 12 system composed of 5 
coupled chains with alternating U = —3, 8, —3, 8, —3. As before, the density of particles 
on the U = 8 chains was kept close to (n) = 1, making them insulating, while the 
density of particles on the U = —3 (superconducting) chains was varied. As before, the 
hopping parameters are t = t' = 1. A pair field A = 0.1 was applied on the bottom 
superconducting chain, and the induced superconducting order parameter was measured 
across the system. Up to m = 2300 states were kept. Fig. [3] shows the induced pair 
fields and the expectation value of throughout the system in two simulations, in 
which the average density of particles on the superconducting chains was (n^c) = 0.7, 
0.47. As expected according to the phase diagram in Fig. [2], in the (usc) = 0.7 case the 
order parameter changes sign from one superconducting chain to the next, while in the 
(^sc) = 0.47 case the sign is uniform. It therefore seems very likely that under the right 
conditions, the two dimensional alternating chain model forms a striped superconductor. 
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Figure 3. The average pair field (/i ~ (c|c|) and spin (5'^) measured in DMRG 
calculations for the 5 x 12 systems with alternating as described in the text. The 
size of the circles indicate the magnitude of (/i, and their color indicate its sign (bright- 
positive sign, dark-negative sign.) The arrows indicate the magnitude and sign of (5^). 
In each calculation, a positive pair field A = 0.1 was applied to the lower chain, and 
a Zeeman field /i = 0.1 was applied to the leftmost site of the second row from the 
bottom. 



5.3. Quasiparticle spectrum of a striped superconductor 

The quasiparticle spectrum of a uniform superconductor is typically either fully gapped, 
or gapless only on isolated nodal points (or nodal lines in 3D). This is a consequence of 
the fact that, due to time reversal symmetry, the points k and — k have the same energy. 
Since the order parameter carries zero momentum, any point on the Fermi surface is 
thus perfectly nested with its time reversed counterpart, and is gapped unless the gap 
function Aj^ vanishes at that point. 

For a striped superconductor, the situation is different. Since the order parameter 
has non-zero momentum Q, only points that satisfy the nesting condition ej^ = e.k+q, 
where ej^ is the single particle energy, are gapped for an infinitesimally weak order. 
Therefore, generically there are portions of the Fermi surface that remain gapless |78j . 
This is similar to the case of a CDW or SDW, which generically leave parts of the 
(reconstructed) Fermi surface gapless, until the magnitude of the order parameter 
reaches a certain critical value. The spectral properties of a striped superconductor 
where studied in detail in Refs. [791 180] . 

As an illustration, we present in Fig. H] the spectral function A(k, = 0) of a 
superconductor with band parameters fitted to the ARPES spectrum of LSCO [45j and 
a striped superconducting order parameter with a single wavevector Q = (27r/8,0) of 
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Figure 4. (Color online.) (a) The spectral function A(k, w = 0) for a striped 
superconductor. The band parameters used in the calculation where fitted to the 
ARPES spectrum of LSCO [45]: t = 0.25, t' = -0.031863, t" = 0.016487, 
t'" — 0.0076112, where t, t',... are nearest neighbor hopping, second-nearest neighbor 
hopping and so on, chemical potential = —0.16235. (All the parameters above are 
measured in eV.) The striped superconducting order parameter has a wavevector of 
Q = (27r/8, 0), and its magnitude is Aq = 60meV. The order parameter is of "d-wave 
character" , in the sense that it is of opposite sign on x and y oriented bonds. The 
thin solid line shows the underlying bare Fermi surface, and the dotted line shows the 
Fermi surface in the presence of the PDW. (b) ^(k, uj) for the same model parameters 
along a cut in k-space. 



magnitude Aq = 60meV. Note that a portion of the Fermi surface around the nodal 
(diagonal) direction remains ungapped (a "Fermi arc" [811 [82]), while both antinodal 
directions [around (vr, 0) and (0, tt)] are gapped. The Fermi arc is in fact the back side of 
a reconstructed Fermi pocket, but only the back side has a sizable spectral weight |83j . 
Its length depends on the magnitude of the order parameter: the larger the magnitude 
of Aq, the smaller is the arc. Note that A(k, u; = 0) is not symmetric under rotation by 
7r/2, because the striped superconducting order breaks rotational symmetry. However, in 
a system with an LTT symmetry (such as LBCO near a; = 1/8 doping) both orientations 
of stripes are present, and an ARPES experiment would see the average of the picture 
in Fig. (Hand its rotation by tt/2. 

6. Order parameter theory of the PDW state 

In this section, we explore the aspects of the theory of a PDW that can be analyzed 
without reference to microscopic mechanisms. We focus on the properties of ordered 
states at T = 0, far from the point of any quantum phase transition, where for the 
most part fluctuation effects can be neglected. (The one exception to the general rule 
is that, where we discuss effects of disorder, we will encounter various spin-glass related 
phases where fluctuation effects, even at T = 0, can qualitatively alter the phases.) For 
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simplicity, most of our discussion is couched in terms of a Landau theory, in which the 
effective free energy is expanded in powers of the order parameters; this is formally not 
justified deep in an ordered phase, but it is a convenient way to exhibit the consequences 
of the order parameter symmetries. 

6.1. Order parameters and symmetries 

We will now define the various order parameters introduced in this section and discuss 
their symmetry properties. The striped superconducting order parameter Aq is a charge 
2e complex scalar field, carrying momentum Q. To define it microscopically, we write 
the superconducting order parameter as 



where R = (r + r') /2, F (r — r') is some short range function (for a "d-wave-like" 
striped superconductor, F (r) changes sign under 90° rotation), and Aq is the uniform 
Q = component of the order parameter. f In the rest of this subsection, we set Aq = 0. 
The effect of Aq is discussed in subsection 16.31 To be concrete, we assume that the 
host crystal is tetragonal, and that there are therefore two potential symmetry related 
ordering wave vectors, Q and Q, which are mutually orthogonal, so Aq must be treated 
on an equal footing with Aq. (The discussion is easily generahzed to crystals with other 
point-group symmetries.) Similarly, for simplicity, spin-orbit coupling is assumed to be 
negligible. 

The order parameters that may couple to Aq and their symmetry properties are 
as follows: The nematic order parameter is a real pseudo-scalar field; the CDW 
Pk with K = 2Q is a scalar field; Sq is a neutral spin-vector field. All these order 
parameters are electrically neutral. Under spatial rotation by tt/2, N — > — A^, pk Pk; 
Sq —>■ Sq, and Aq ±Aq, where ± refers to a d-wave or s-wave version of the striped 
superconductor. Under spatial translation by r, A^ — > A^, px ^^^ ^Pk, Sq c'^^ ^Sq, 
and Aq c'^ '^Aq. Note that since the SDW and CDW orders are real, Sq = S-q 
and = p-K- Generally, Aq and Aq are independent. 

f A state in which both components of the SC order parameter coexist, Aq 7^ and Aq ^ is certainly 
not "uniform" . Even a weak Aq ^ imphes the existence of a modulation of the local amplitude 
of the SC order parameter, and a SC state is "truly uniform" only if Aq — 0. Nevertheless, as we 
will see in Section [6.31 the properties of a SC state in which both order parameters coexist are largely 
dominated by the "uniform" component Aq, and the striking features of the PDW state are not directly 
observable. In this sense, the uniform-PDW coexisting SC state is effectively "uniform." 
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6.2. Landau Theory 

Specifically, the emphasis in this section is on the interrelation between striped 
superconducting order and other orders. There is a necessary relation between this 
order and CDW and nematic (or orthorhombic) order, since the striped superconductor 
breaks both translational and rotational symmetries of the crystal. From the microscopic 
considerations, above, and from the phenomenology of the cuprates, we also are 
interested in the relation of superconducting and SDW order. The Landau effective 
free energy density can then be expanded in powers of these fields: 

= J^2 + ^3 + -5^4 + ... (19) 

where the quadratic term, is simply a sum of decoupled terms for each order 
parameter, 

-^3 = 1s[P-kSq ■ Sq + p_kSq ■ Sq + c.c] (20) 

+ 7a[P-kA^_qAq + p_kA^qAq + c.c] 

+ (7AiV[A^AQ + A*_qA Q - A^Aq - A* qA q] 

+ gsN[S^Q ■ Sq - S Q ■ Sq] 

+ 9cN[p.KpK - P-kPk], 
and the fourth order term, which is more or less standard, is shown explicitly below. 

The effect of the cubic term proportional to 7^ on the interplay between the spin and 
charge components of stripe order has been analyzed in depth in [HI]. Similar analysis 
can be applied to the other terms. In particular, the 7a and terms imply| that the 
existence of superconducting stripe order (Aq 7^ 0, and Aq = 0), implies the existence 
of nematic order (A^ 7^ 0) and charge stripe order with half the period {p2Q 7^ 0). 
However, the converse statement is not true: while CDW order with ordering wave 
vector 2Q or nematic order tend to promote PDW order, depending on the magnitude 
of the quadratic term in JF2, PDW order may or may not occur. 

One new feature of the coupling between the PDW and CDW order is that it 
produces a sensitivity to disorder which is not normally a feature of the superconducting 
state. In the presence of quenched disorder, there is always some amount of spatial 
variation of the charge density, p(r), of which the important portion for our purposes 
can be thought of as being a pinned CDW, that is, a CDW with a phase which is a 
pinned, slowly varying function of position, p(r) = |pk| cos[K ■ r + (/>(r)]. Below the 
nominal striped superconducting ordering temperature, we can similarly express the 
PDW order in terms of a slowly varying superconducting phase, A(r) = |Aq| exp[zQ • 
^ + ^^Q(r)] + |A_q| exp[— iQ ■ r + i6'_Q(r)]. The resulting contribution to JF3 is 

^3,7 = 27a|pkAqA_q| cos[2^^_(r) - 0(r)]. (21) 

I Note that the 7a term is odd under a particle-hole transformation, which takes pK — > — Pk- 
Therefore, if the system has exact particle-hole symmetry, this term vanishes, and there is no necessary 
relation between Aq and pK- Microscopic systems are generically not symmetric under charge 
conjugation. However, some real systems {e.g. the cuprates) are not too far from being particle-hole 
symmetric, and therefore in these systems 7a is expected to be relatively small. 
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where 

0±ir) =[^^Q(r)±^^_Q(r)]/2; (22) 
e±Q(r) = [e+(r)±0_(r)]. 

The aspect of this equation that is notable is that the disorder couples directly to a 
piece of the superconducting phase, 9-. No such coupling occurs in usual momentum 
superconductors. 

It is important to note that the condition that A(r) be single valued implies that 
^Q(r) and ^^-Q(r) are defined modulo 2tt. Correspondingly, 9± are defined modulo vr, 
subject to the constraint that if 6± —>■ 6± + TT'm± then m+ + m_ must be an even integer. 
Since (p and 9_ are locked to each other at long distances, the possible topological 
excitations of the coupled PDW-CDW system are thus point defects in 2D and line 
defects in 3D classified by the circulation of 6+ and on any enclosing contour. The 
elementary topological defects thus are: a) An ordinary superconducting vortex, about 
which = 2TT and A(f) = 0. b) A bound-state of a half vortex and a dislocation, § 
about which A^^+ = vr and A(f) = 27r. c) A double dislocation (or dislocation bound 
state) about which A6+ = and A0 = Att. All these defects have a logarithmically 
divergent energy in 2D, or energy per unit length in 3D; the prefactor of the logarithm 
is determined by the superfiuid stiffness for the vortex, the elastic modulus of the CDW 
for the double vortex, and an appropriate sum of these two stiffnesses for the half vortex. 
Consequences of this rich variety of topological defects are discussed in |6l [25|, [86] 

An important consequence of the coupling between the superconducting and CDW 
phase is that the effect of quenched disorder, as in the case of the CDW itself, destroys 
long-range superconducting stripe order. (This statement is true [87], even for weak 
disorder, in dimensions d < A.) Naturally, the way in which this plays out depends on 
the way in which the CDW state is disordered. 

In the most straightforward case, the CDW order is punctuated by random, pinned 
dislocations, i. e. 2ti vortices of the field. The existence of the coupling in Eq. [2T] implies 
that there must be an accompanying vr vortex in 9^. The condition of single- valued-ness 
implies that there must also be an associated half- vortex or anti- vortex in 6'+. If these 
latter vortices are fluctuating, they destroy the superconducting state entirely, leading 
to a resistive state with short-ranged striped superconducting correlations. If they are 
frozen, the resulting state is analogous to the ordered phase of an XY spin-glass: such 
a state has a non-vanishing Edwards-Anderson order parameter, spontaneously breaks 
time-reversal symmetry, and, presumably, has vanishing resistance but no Meissner 
effect and a vanishing critical current. In 2D, according to conventional wisdom, a spin- 
glass phase can only occur at T = 0, but in 3D there can be a finite temperature glass 
transition |88j . 

In 3D there is also the exotic possibility that, for weak enough quenched disorder, 
the CDW forms a Bragg-glass phase, in which long-range order is destroyed, but no free 

§ The possibility of half vortices in a striped superconductor and their effect on the phase diagram in 
the clean case was discussed by D. F. Agterberg and H. Tsunetsugu[55]. 
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dislocations occur [211 [101 EI]- In this case, can be treated as a random, but single- 
valued function - correspondingly, so is 9_. The result is a superconducting Bragg-glass 
phase which preserves time reversal symmetry and, presumably, acts more or less the 
same as a usual superconducting phase. It is believed that a Bragg-glass phase is not 
possible in 2D [90]. 

Another perspective on the nature of the superconducting state can be obtained 
by considering a composite order parameter which is proportional to AqA_q. There is 
a cubic term which couples a uniform, charge 4e superconducting order parameter, A4, 
to the PDW order: 

J^;, = (74{A*[AqA_q + AqA_q] + c.c.} (23) 

This term implies that whenever there is PDW order, there is also necessarily charge 
4e uniform superconducting order. However, since this term is independent of it 
would be totally unaffected by Bragg-glass formation of the CDW. The half-vortices 
in 6'_|_ discussed above can simply be viewed as the fundamental (hc/Ae) vortices of a 
charge 4e superconductor. 

Some additional physical insight can be gained by examining the quartic terms (^4 
in Eq. [T9l) . Let us write all the possible fourth order terms consistent with symmetry: 

^4 = M (S^Q ■ S'qA^A^q + ■ S'qA^A.q + c.c.) 

+ (v+[S^Q ■ Sq + S^Q ■ Sq] + v+[\pk\^ + IpkH) 
X (|AqP + |A_qP + |AqP + |A_qP} 

+ ■ 5q - Sa ■ Sa] + V-[\Pk\' - \Pk\']) 

x(|AQp + |A_Qp-|AQp-|A_Qn 
+ .AT^ {(|AqP + IA.qH + (|AqP + IA^qH} 

+ { (|AqP + |A_Qn V (IAqI^ + IA^qH^} 

+ A- {(|AqP - |A_Qn% (|AqP - IA^qH^} 
+ A(|AqP + |A_qP)(|Aq|2 + |A_qP) 

+ . . . (24) 

where we have explicitly shown all the terms involving Aq, while the terms . . . represent 
the remaining quartic terms all of which, with the exception of those involving N, are 
exhibited explicitly in [8^ . 

There are a number of features of the ordered phases which depend qualitatively on 
the sign of various couplings. Again, this is very similar to what happens in the case of 
CDW order - see, for example, [921 [93]. For instance, depending on the sign of A, either 
unidirectional (superconducting stripe) or bidirectional (superconducting checkerboard) 
order is favored. 

On physical grounds, we have some information concerning the sign of various 
terms in JF4. The term proportional to u determines the relative phase of the spin 
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and superconducting stripe order — we believe u > which thus favors a 7r/2 phase 
shift between the SDW and the striped superconducting order, i.e. the peak of the 
superconducting order occurs where the spin order passes through zero. The other 
interesting thing about this term is that it imphes an effective cooperativity between 
spin and striped superconducting order. The net effect, i.e. whether spin and striped 
superconducting order cooperate or fight, is determined by the sign of |u| — f + — f _, such 
that they cooperate if \u\ — — v_ > and oppose each other if |u| — 1>+ — < 0. It is 
an interesting possibihty that spin order and superconducting stripe order can actually 
favor each other even with all "repulsive" interactions. The term proportional to A_ 
determines whether the superconducting stripe order tends to be real (A_ > 0), with a 
superconducting order that simply changes sign as a function of position, or a complex 
spiral, which supports ground-state currents (A_ < 0). 

6. 3. Coexisting uniform and striped order parameters 

Finally, we comment on the case of coexisting uniform and striped superconducting 
order parameters. Such a state is not thermodynamically distinct from a regular 
(uniform) superconductor coexisting with a charge density wave, even if the uniform 
superconducting component is in fact weaker than the striped component. Therefore, we 
expect many of the special features of the striped superconductor (such as its sensitivity 
to potential disorder) to be lost. Here, we extend the Landau free energy to include a 
uniform superconducting component, and show that this is indeed the case. 

We will now analyze the coupling of a striped superconducting order parameter 
Aq to a uniform order parameter, Aq. In this case, we have to consider in addition to 
the order parameters introduced in Sec. [6] a CDW order parameter with wavevector Q, 
denoted by pq. The additional cubic terms in the Ginzburg-Landau free energy are 

J^3,„ = 7qAo [pqA_q + p_qAq + PqA_q + P_qAq] + c.c. 



+ 9p 



(25) 



P-2qPq + P-2qPq + c.c. 

Eq. [25] shows that if both Aq and Aq are non-zero, there is necessarily a coexisting 
non-zero pq, through the 7q term. The additional quartic terms involving Aq are 



^4,« = ^^A (A^2AqA_q + A^2AqA_q + c.c.) + (5|Ao|'[|AqP + |A 



Ql 



1 21 



+ |Ao| 



Up (\pQ? + ipqT) + (ip2qi^ + iP2Qr 



+ i;'|Aon5_Q-5Q + 5_Q-^Q]. (26) 

Let us now consider the effect of quenched disorder. Following the discussion preceding 
Eq. [221 we write the order parameters in real space as 

A (r) = |Ao| e'^« + |Aq| e^'^'^+Q") + |A_q| e^'^-^^"^") (27) 

and 

P{f) = Ipq I cos (Q ■ r + (/)q) + |p2Q I cos (2Q ■ r + (f)) . (28) 
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Let us assume that the disorder nucleates a point defect in the CDW, which in this 
case corresponds to a 27r vortex in the phase </)q. By the Qp term in Eq. [251 
this induces a Att vortex in 0. (Note that in the presence of pq, a 2tt vortex in 
is not possible.) The 7a term in Eq. [21] then dictates a 2it vortex in the phase 
6'_ = (^^Q — 9-q) /2. However, unlike before, this vortex does not couple to the 
global superconducting phase 9+ = {9q + 9_q)/2. Therefore, an arbitrarily small 
uniform superconducting component is sufficient to remove the sensitivity of a striped 
superconductor to disorder, and the system is expected to behave more or less like 
a regular (uniform) superconductor, albeit with a modulated amplitude of the order 
parameter. 

Since the usual (uniform) superconducting order and the PDW break distinct 
symmetries, nothing can be said, in general, about the conditions in which they 
will coexist. However, microscopic considerations can, in some cases, yield generic 
statements, too. For example, in a striped SC, a uniform component of the order 
parameter can be generated by dimerizing the stripe order, such that the positive and 
negative strips of superconducting order are made alternately broader and narrower. In 
any structure (such as the LTT structure of LBCO), in which there is zero Josephson 
coupling between neighboring layers, a coupling is generated, thus lowering the energy 
of the system, in proportion to the square of the dimerization. Presumably, so long 
as the PDW period is incommensurate with the underling lattice, there is also a 
quadratic energy cost to dimerization which is related to an appropriate generalized 
elastic constant of the PDW. However, if the PDW has a long period, this elastic 
constant will be vanishingly small. Thus, any long period, incommensurate PDW may 
generically be expected to be unstable toward the generation of a small amount of 
uniform SC order. 

7. Non-collinear order and time reversal symmetry breaking 

In a layered system, PDW order in the planes can lead to frustration of the inter-plane 
Josephson coupling, which naturally explains the layer decoupling seen in 1/8 doped 
LBCO.f In analogy with frustrated magnetic systems (in which the superconducting 
order is thought of as an XY pseudo-spin), this frustration can also lead to various 
forms of non-collinear order. In the PDW case, such non-collinear orders break time- 
reversal symmetry and are accompanied by spontaneous equilibrium currents. 

In this section, we give detailed predictions for the patterns of bulk time-reversal 
symmetry breaking and spontaneous currents in various lattice geometries. We will 
discuss this problem at zero temperature and at a classical level. It is worth noting that 
the non-collinear order, where it occurs, results in a partial lifting of the frustration. 
In the case of a PDW in the LTT structure (relevant to La2_a;Ba2:Cu04), we shall 
show that it results in a non-vanishing effective Josephson coupling between planes, 

f The problem of the 3D phase transition in a system with an effective layer decoupling is largely 
unsolved. See, however, the recent work of Raman, Oganesyan and Sondhi.[94] 
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Figure 5. (Color online.) (a) Model for a striped superconductor with an LTT 
structure. Solid (dashed) lines represent positive (negative) Josephson couplings. The 
arrow on the center of each link indicates the direction of the equilibrium current across 
that link. The red arrows on the vertices represent the superconducting phases, (b) 
Same as (a) for an orthorhombic striped superconductor, where the charge stripes are 
shifted by half a period from one layer to the next, (c) An in-plane domain wall. 

and hence, in a sense, spoils the strict layer decoupling we have touted. However, this 
effective Josephson coupling is equivalent to a higher order coupling [6j (due to coherent 
tunneling of two Cooper pairs), both in terms of its small magnitude, and its dependence 
on the cosine of twice the difference of the superconducting phases on neighboring planes. 
(See Eqs. [30] and [321) Note also that defects (such as point defects, domain walls or 
twin boundaries) can lead to additional intra-plane time reversal symmetry breaking, 
that can drive the system into a glassy superconducting state (as discussed in Sec. [6]).$ 
Let us start with the case of the LBCO LTT structure, in which the stripe direction 
rotates by 90° between adjacent planes. We model the system by a 3D discrete lattice 
of Josephson junctions, shown in Fig. [5tL.§ The lattice spacing in the plane is the 
inter-stripe distance A, and c is the inter-plane distance. Each lattice point has a single 
degree of freedom which is the local value of the superconducting phase at that point. 

J An in-plane magnetic field can also change the inter-layer frustration, leading to small violations of 
the layer decoupling effect. If large enough such effects can be used to detect a PDW state. A similar 
effect can also take place in junctions between an FFLO state and a uniform superconductor |95j . 
§ Note that we are actually considering a simplified version of the LBCO LTT structure. The structure 
in Fig. [SK has two planes per unit cell, while the LBCO LTT structure has four. The difference is that 
in LBCO, the charge stripes in second neighboring planes (which are parallel to each other) are shifted 
by half a period relative to one another, while in Fig. [5^ they are not. However, the considerations we 
discuss here are the same for two structures, and the resulting non-coUinear ground states are similar. 
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J, — J', J" are the intra-stripe, the inter-stripe and the inter-plane Josephson couphngs, 
respectively. We assume that J > J' ^ J" > 0, corresponding to a unidirectional 
striped superconductor in the planes. For any collinear configuration, the Josephson 
coupling between the planes vanishes. However, if the staggered order parameter in 
each plane is rotated by 90° relative to its neighbors, then the energy can be lowered 
by distorting the phases periodically with respect to the collinear configuration in each 
plane. We use a variational ansatz for the phases 6^ of the form 

= ^L±tlLyT, + l^i^ (^ + 1)^ + (29) 



2 2 ^ ^ 

where r = (x, y, z) is the integer valued position vector (x and y are measured in units of 
A, and z is measured in units of c), and the distortion angle 6' is a variational parameter. 
The Josephson energy per site function of 6 is 

Eltt {0) = -{J+ f) cos 26 - J" sin 29. (30) 

The inter-plane coupling energy gain is linear in 9, whereas the cost in intra-plane 
coupling energy is quadratic in 9. Thus the distortion occurs for any non-zero value of 
the inter-plane coupling J" . Minimizing Eq. (!30|) . we get 

/" 

tan2^ = - -. (31) 

J + J' ^ ' 

The equilibrium currents across the three types of links are J = J sin 20, J'' = J' sin 29 

and J'" = J" cos 29 = J' + J'', where Eq. l3T]was used in the last relation. The directions 

of the currents are as indicated in Fig. [5^. Associated with these currents is a magnetic 

field with non-zero components in all three directions. The wavevector associated with 

this pattern is Q = j, where A is the inter-stripe distance (for LBCO at x = 1/8, 

A ~ 4a where a is the Cu-Cu distance) and c is the inter-plane distance. 

The non-coUinear distortion in the above pattern induces an effective non-zero 

inter-plane coupling. However, the effective inter-layer coupling is (taking the limit 

J" < J, J') 

and is therefore much smaller than the bare inter-plane coupling J". Note, moreover, 
that the induced Josephson coupling between two neighboring planes with PDW 
superconducting phases 9i and 9j has the form Jeg cos[2(6'j — 9j)], i.e. its period in 
the relative phase is tt. 

Next, we consider the case of an orthorhombic structure (such as the LTO phase 
of LBCO). In this case, rotational symmetry in the plane is broken in the same way in 
every plane, and the stripes are all in the same direction. However, we assume that due 
to Coulomb interactions, the charge stripes are shifted by half a period between adjacent 
planes. (Such a shift is indeed observed between second neighbor planes in the LTT 
phase of LBCO, in which the stripe direction is parallel.) Therefore, the inter-plane 
coupling is frustrated due to the resulting "zigzag" geometry. We shall show below that 
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the ground state has spiral order which partially relieves this frustration. Introducing a 
spiral twist angle 9, such that 9j. = 2x9 (as shown in Fig. [Sjo), costs an energy -Eqrt (9) 
per stripe, given by 

EoRT {9) 



L 



J' cos 29 -2 J" cos 9 (33) 



where L is the length of each stripe. The minimum is for cos 6* = Therefore 
a spiral distortion occurs for any J" < 2J'. The currents along this links are 
J'" = —J' = J" sin 9, and their directions are indicated in Fig. [Sb- Each plane 
carries a uniform current which flows perpendicular to the stripes, and an equal and 
opposite current flows between the planes. The magnetic field associated with these 
currents is pointing parallel to the stripe direction, and its lowest Fourier component is 
at wavevector Q = (O, 0, ^) . 

Finally, we turn to the case of a domain wall in the PDW order, depicted in Fig. 
|5t. (Such a defect is very costly energetically, but it is favored by a twin boundary 
in the crystal structure.) The Josephson coupling across the domain wall vanishes 
for any collinear configuration. The energy can be lowered by distorting the phases 
in the pattern shown in Fig. [5b, which is closely analogous to the minimum energy 
configuration in the LTT case (Fig. [5K)- The superconducting phases 9r are given by 

_ J (a; + i)vr-(-l)^^. (x<l) 
^ " \ yvr + {-l)y9. (x > 1) ' 

where the distortion angle 9^ depends on the distance from the domain wall x. (In our 
notation, x = and 1 are the two columns on either side of the domain wall.) The 
energy is 

^I^l^M^ = - J cos {9.^, - 9^) -J'f2 cos (29.) 

x=l x=l 


-J' Yl cos{9,^i-9^)-J J2 cos (20,.) 

x=—oo x=—oo 

- Jsin(0,=i + 0,=o). (35) 

Here, J is the Josephson coupling across the domain wall, and L is the number of sites 
along the domain wall. For simplicity, we assume that J <^ J, J', in which case 9x ^ 
and we may expand Eq. (1351) to second order in 0,. Minimizing -Edis ({^a;})) obtain 
the following solution: 

9 - I ^<""' r3fii 

- ^ 0>e-/^^ (x > 1) ' ^^^^ 



where a = 2sinh-^f./4\ /? = 2 sinh"^ f ./^'l , 0< = and 9 



J 



^+47- ■ Associated with the distortion of the superconducting phases is a periodic 
pattern of spontaneous currents, shown in Fig. [5t, with periodicity of two inter-stripe 
distances. 
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Similar considerations apply to an in-plane Josephson junction between a striped 
superconductor and a uniform superconductor, if the boundary is perpendicular to the 
stripe direction. Therefore, in such a junction time reversal symmetry is also broken. 
The critical current is of order . . [This follows from the same considerations as 

mm|J,J'| L 

the effective inter-plane coupling in the LTT case, Eq. [321] It is thus suppressed relative 
to the critical current of a Josephson junction between uniform superconductors, which 
is of order J, as a result of the frustration of the Josephson coupling across the junction. 
Similarly to inter-plane coupling in the LTT case, the period of the coupling between a 
uniform and a striped superconductor in the relative phase is vr, i.e. half of the period 
of the coupling between two uniform superconductors. 

8. Connections and History 

The notion of a superconducting state with spontaneously generated oscillations in 
the sign of the order parameter has cropped up, under various guises, a number of 
times in the past. It is worthwhile to recount some of these circumstances, not only in 
the interest of scholarship, but also to broaden the range of phenomena which can be 
addressed within the same conceptual framework. 

8.1. Josephson IT junctions 

Since the superconducting order parameter is a charge 2e scalar field, it is often assumed 
that it is possible to think of the superconducting state as a Bose condensed state of 
charge 2e bosons. In contrast, most classic treatments of the subject |96] emphasize that 
many features of BCS theory, especially those associated with quasiparticle coherence 
factors, cannot be understood in this way. At the very least, a bosonic theory is 
inadequate to capture basic features of the groundstate of any superconductor which has 
gapless quasiparticles, either because of the order parameter symmetry {e.g. d-wave) or 
because of scattering from magnetic impurities (gapless superconductor). 

Even ignoring the possibihty of gapless quasiparticles, there are qualitative 
possibilities in a fermionic system that cannot occur in a bosonic system. A feature 
of a time reversal invariant bosonic system is that the ground-state can be chosen to 
be real and nodeless. Thus, the order parameter in a Bose-condensed system must 
have a phase which is independent of position. The tt junctions, which we have been 
discussing, are possible only because of the composite character of the superconducting 
order parameter [71] . 

There have been several previous theoretical studies which have found 
circumstances under which vr junctions might occur [7T1 [72], [97]. More recently, the 
existence of such vr junctions in the predicted circumstances have been confirmed 
by experiment. The first such experiments [711 ES] were significant as the "phase 
sensitive" measurements which definitively established the d-wave symmetry of the 
superconducting order in the cuprates. More recently, however, mesoscopic vr 
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junctions between two s-wave superconductors have been constructed and characterized 
[73] . In our opinion, these latter experiments are also landmarks in the study of 
superconductivity. They establish that vr junctions, the essential ingredient for the 
existence of striped superconductors, are physically possible. 

8.2. FFLO states 

In a superconductor with negligible spin-orbit coupling, it is possible to generate an 
imbalance in the population of up and down spin quasiparticles, either by applying 
a magnetic field in a geometry in which it predominantly couples to the electron 
spins, or by injecting a non-equilibrium population of quasiparticles from a neighboring 
ferromagnet [9Hj. In the related systems of cold fermionic atomic gases, it is possible 
to vary the population of up and down spin atoms independently, and to study the 
effect of this population imbalance on the superfluid state [99], llOOl 125] . While a first 
order quenching of the superconducting state is possible under these circumstances, 
there has also been considerable discussion of the possibihty of spatially modulated 
superconducting states, so called FFLO states [TT|,[T2]. Two distinct states of this sort 
have been considered: (1) The FF state [11], in which the order parameter has constant 
magnitude but a phase which twists as a function of position according to ^ = Akp ■ r, 
where Ahp is the difference between the up spin and down spin Fermi momentum. (2) 
The LO state [12], in which the order parameter remains real, but oscillates in sign with 
a period L = 27r/|Aki?|. 

The LO state is similar in structure to the striped superconductor considered here. 
In the order parameter theory presented in Sec. 16.21 it corresponds to A_ > in Eq. 
(!24|) . The parallel with the FF state (which is realized in the order parameter theory 
for A_ < 0) is less crisp, but when superconducting striped spirals which spontaneously 
break time reversal symmetry arise due to the appropriate type of geometric frustration 
of the Josephson couplings (as discussed in Sec. [7]), states that are in many ways 
analogous to the FF state also occur in striped superconductors. Thus, many of the 
physical phenomena we have discussed in this paper are pertinent to the FFLO phases in 
more weakly correlated systems, with the added richness p5] in the case of cold atomic 
gasses that there are conserved quantities associated with the continuous rotational 
invariance of the underlying Hamiltonian. 

However, the FFLO states arise from the explicit breaking of time reversal 
symmetry. Absent a magnetic field, Kramer's theorem implies perfect nesting between 
time-reversed pairs of states on opposite sides of the Fermi surface, so BCS pairing 
always occurs preferentially at k = 0. This constraint is removed when time reversal 
symmetry is explicitly broken. One can think of the FFLO states as taking advantage 
of the "best" remaining approximate nesting vector, Akp, in the two-particle channel. 
Alternatively, one can think of the LO state as consisting of a set of discommensurations 
[981 [25] such that the excess spin-up quasiparticles are incorporated in mid-gap states 
localized near the core of the discommensuration. 
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The energetic considerations that lead to the FFLO states are thus very different 
than the strong-couphng physics that gives rise to the striped superconductor. f The 
fact that the FFLO states exphcitly break time reversal symmetry implies that they are 
macroscopically distinct (as phases of matter) from the striped superconductors that 
preserve this symmetry. Even in comparison with striped states which spontaneously 
break time reversal symmetry, the distinction remains that the FFLO states have a 
net magnetization, while the striped superconductor does not. Conversely, the FFLO 
states generally have no particular relation to other flavors of electronic ordering, while 
striped superconductors, as is characteristic of all electronic liquid crystals, embody a 
subtle interplay between multiple ordering tendencies. Specifically, since the striped 
superconductor seems to be generally associated with the strong coupling physics of 
doped antiferromagnets, there is a natural sense in which antiferromagnetism, charge 
density wave formation, and striped superconductivity are intertwined. 

8.3. Intertwined orders and emergent symmetries 

One explicit way in which the relation between several order parameters can be more 
intimate than in a generic theory of "competing orders" is if there is an emergent 
symmetry at low energies which unifies them. In particular, the order parameter 
structure of the PDW state, involving several order parameters coupled to each 
other, evokes the S0{5) approach of a unified description of antiferromagnetism and 
uniform c?-wave superconductivity flSl E]. Indeed, by tuning the parameters of the 
effective Landau-Ginzburg theory that we presented in other sections it is possible to 
achieve an effective enlarged symmetry which makes it possible "rotate" the striped 
superconducting order and charge stripe order parameters into each other. Even if 
the enlarged symmetry is not exact, a rotation of the order parameters is possible 
but with a finite energy cost ( similar to a "spin flop".) It is also worth noting that 
a symmetry which allows a similar form of unification of d-wave superconductivity, 
electron nematicity, and (i-density wave order [61] (dDW) has recently been found to 
exist under special circumstances by Kee et al. |103j . It is therefore possible that there 
could exist additional forms of striped superconducting states which interleave these 
orders. 

Thus, it is possible to view the PDW state as a "liquid crystalline" analog of the 
5*0(5) scenario. Indeed, the possibility of an 50(5) "spiral" was discussed previously 
by Zhang [104j . However, it should be noted that in the context of any conventional 
Landau-Ginzburg treatment of a system of competing orders, a general theorem [105] 
precludes a sign change of any component of the order parameter, and hence precludes 
the existence of spirals. In order to get a PDW state from an interplay between d- 
wave superconductivity and antiferromagnetism, unconventional gradient dependent 
interactions between the different order parameters, such as those discussed in |105j . 



f FFLO states in the absence of magnetic fields have been shown to exist for special band structures 
in lD[ini] and 2D|Tg^. 
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must play a significant role in the physics. 

In other words, in addition to the standard couplings allowed by a theory with 
several order parameters, the existence of a stripe order (for instance) in the charge 
order parameter must be able to induce a texture in the superconducting order as well. 
A useful analogy to keep in mind is the McMillan-deGennes theory of the nematic- 
smectic transition in classical liquid crystals in which the nematic order parameter acts 
as a component of a gauge field thus coupling to the phase of the smectic order, or in 
blue phases of liquid crystals. (For a detailed discussion of these topics in liquid crystals 
see, e.g. |106l 1107] .) In fact, Ref [25] presents a theory of FFLO states in ultra-cold 
atoms with gauge-like couplings [i.e. covariant derivative couplings) that relate the 
stripe (and spiral) order to the superconducting order. 

In addition to the conceptual advantages, noted above, the liquid-crystal picture of 
the PDW state offers a direct way to classify the phase transitions (both quantum and 
thermal) out of this state. Thus, in addition to a direct transition to a normal state, 
intermediate phases characterized with composite order parameters, are also possible 
leading to an interesting phase diagram. We will explore these issues in a separate 
pubhcation [86] . 

8.4. PDW states in Hubbard and t-J models 

In the context of the cuprates, there have been several studies looking for a striped 
superconducting state in the t — J or Hubbard models. On the one hand, extensive, but 
not exhaustive DMRG calculations by White and Scalapino [771 1108] have consistently 
failed to find evidence in support of any sort of spontaneously occurring vr junctions. 
On the other hand, a number of variational Monte Carlo and renormalized mean field 
calculations have concluded that the striped superconductor is either the ground-state 
of such a model [109] . under appropriate circumstances, or at least close in energy to 
the true ground state |110lllll[ [80]. These latter calculations are certainly encouraging, 
in the sense that they suggest that there is no obvious energetic reason to rule out 
the existence of spontaneously occurring PDW order in strongly correlated electronic 
systems. However, the fact remains that no spontaneous vr-junction formation has yet 
been observed in DMRG or other "unbiased" studies of the t — J ot the purely repulsive 
Hubbard models, indicating that there remain basic unsettled issues concerning the 
microscopic origins of vr junctions. 

9. Final thoughts 

In this paper, we have introduced the PDW phase and studied its properties 
theoretically. In terms of symmetry, the PDW is distinct from the standard uniform 
superconductor. While some of its properties are similar to those of a uniform 
superconductor {e.g., zero resistance), others are markedly different: most importantly, 
the existence of a Fermi surface (and hence a finite density of states) in the ordered phase 
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[inilSS], the possibility of frustration of the inter-layer coupling (depending on the lattice 
geometry), and the strong sensitivity to (non-magnetic) disorder. Generically, the PDW 
state in the presence of weak disorder is expected to give way to a "superconducting 
glass" phase, in which the configuration average of the local superconducting order 
parameter vanishes, but the Edwards- Anderson order parameter is non-zero (and hence 
gauge symmetry is broken). 

Even though the ordered PDW state itself is time reversal invariant, time reversal 
symmetry breaking is a very natural consequence of PDW order, either in the 
superconducting glass phase, or as a way of relieving the frustration of the Josephson 
couplings in some crystal structures. Specifically, frustration can lead to non-collinear 
ground state configurations of the superconducting pseudo-spins (representing the local 
phase of the superconducting order), which are analogous to the non-collinear ground 
states which are often found in frustrated spin systems. An even more exotic state that 
can naturally emerge from a "parent" PDW state is a superconductor with a charge 4e 
order parameter [HI [3, [25] , which can result when the CDW part of the PDW order is 
melted by either quantum or thermal fluctuations. 

The occurrence of PDW states in microscopic models is an intrinsically strong 
coupling effect, since PDW order (much like CDW or SDW) is not an instability of 
a generic Fermi surface. In this paper, we have provided a "proof of principle" of a 
not-too-contrived, strongly correlated, microscopic model with a PDW ground state. 
This model mimics some features of the striped state found in the cuprates {e.g., it has 
charge stripes separated by vr-phase-shifted spin stripes). Whether a PDW state can 
be found in more realistic models, which include such features as uniformly repulsive 
interactions and a d-wave-like order parameter, remains to be settled. 

Doped Mott insulators are strongly correlated systems whose ground states have 
a strong tendency to form liquid- crystalline- like |8] inhomogeneous phases. fll2[ I113[ 
11141 11151 I116j . In this regard, the PDW state is an electronic liquid crystal phase 
in which the superconducting and charge/spin orders do not compete with each other 
but rather are intertwined. As some of us have noted earlier, [1171 1118] the observation 
of a high pairing scale in such an electronically inhomogeneous state is suggestive of the 
existence of an optimal degree of inhomogeneity for superconductivity. Indeed, recent 
ARPES data suggest that the stripe order that develops in La2-a:Baa;Cu04 does not 
suppress the pairing scale. f The fact that the pairing scale is large in this material 
suggests that the development of charge stripe order suppresses the development of 
superconducting coherence but not pairing. In fact, it gives credence to the argument 
that there is a connection between the emergence of charge order and the mechanism of 
superconducting pairing. [1171 1118j 

However, at present, it is unclear to what extent PDW order should be expected 
to be common where stripe order occurs. On the purely theoretical side, PDW 
order has proven elusive in DMRG studies of models [77[ llOSj with entirely repulsive 

f ARPES data in La2-a;Baa;Cu04 shows a substantial and weakly doping dependent anti-nodal gap 
accross x = 1/8(111 US], where the signatures of the PDW state are strongest. 
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interactions. Indeed, in a previous publication [7j, we showed that in any weakly 
interacting superconductor, vr junctions can only occur under exceedingly fine-tuned 
circumstances. It is clear from variational calculations [IDHl IllOj lllll 180] that for 
strong interactions, the differences in energy between the PDW and uniform sign 
superconducting states in striped systems is relatively small; what particular features of 
the microscopic physics tip the balance one way or another is still not clear. Accordingly, 
it is not clear, in the absence of unambiguous experimental evidence, whether in the 
context of the cuprates, we should expect the PDW state to be a rare occurrence, perhaps 
stabihzed by some particular detail of the electronic structure of La2_a;Baa;Cu04, or if 
instead we should infer that some degree of local PDW order exists in any cuprate in 
which evidence of local stripe correlations can be adduced. 

To close this Section, we turn to discuss the evidence for PDW states in the cuprate 
high temperature superconductors. The analysis of the PDW state was motivated by the 
experimental observations on La2_a;Baa;Cu04. Having studied the nature of this phase, 
we will now discuss to what extent the signatures of the PDW state are consistent 
with experiment. Finally, we speculate on the possible relevance of these ideas to other 
members of the cuprate family. 

9.1. Striped SC phases in La2-xBaxCu04, and 214 cuprates 

As already discussed in Sec. [3l the onset of clearly identifiable 2D superconducting 
correlations in La2-xBa2:Cu04 with a; = | occurs at ~ 40 K, together with the onset 
of static spin-stripe order. It would be natural to associate this behavior with the 
simultaneous onset of local PDW order; however, an attempt to reach a consistent 
interpretation of a broad range of results leads to a more nuanced story. 

The original motivation for applying the PDW concept to Lai,875Bao.i25Cu04 was 
to explain the dynamical layer decoupling through the frustration of the interlayer 
Josephson coupling in the LTT phase [6], 1109] . as discussed in Sec. [71 It provides a 
compelling account:}: for the induced dynamical layer decoupling produced in underdoped 
La2-xSr2:Cu04 by a modest c-axis magnetic field [2]. Moreover, the sensitivity of 
the PDW to disorder which limits the growth of the superconducting correlation 
length within the planes, provides a natural explanation for the existence of an 
enormously enhanced "superconducting fluctuation" regime, characterized by enhanced 
contributions of local superconductivity to the electrical conductivity and to (strongly 
anisotropic) diamagnetism, but with no global phase coherence. Thus, it naturally 
accounts for the most dramatic aspects of the experimental data [1] below the spin 
ordering temperature Tsdw- We consider this strong evidence that the basic ingredients 
of the theory are applicable to the stripe ordered state of La2_a;Baa;Cu04 and closely 

t Since La2_2:Sr2:Cu04 retains the LTO structure to low temperatures, and the spin correlations in 
the c-direction measured at zero field are extremely short-ranged [119j , it is unclear whether the charge 
stripes in neighboring planes tend to be perpendicular to each other, as in the LTT materials, or 
parallel but ofl^set by half a period from each other, as in the YBa2Cu306+2: bilayers. In either case, 
the interlayer Josephson coupling for a PDW would be highly frustrated. 
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related materials. In addition, the observed transition at temperature T^d into a 
state with zero resistance in all direction has a natural interpretation in terms of 
an assumed PDW state as the superconducting glass transition |7j. Besides having 
zero resistance, the glass phase presumably shows no Meissner effect and zero critical 
current. If this latter identification is correct, it leads to the further prediction that this 
phase should be characterized by various phenomena associated with slow dynamics, 
characteristic of spin glasses, as well as with breaking of time reversal symmetry. The 
experimental detection of such phenomena below T30 supercurrents) would serve as 
further confirmation of the existence of a PDW in this material. (For example, the glass 
phase would likely exhibit a metastable zero- field Kerr effect [58J.) 

One can also look for evidence for the PDW in single-particle properties. One 
of the key features of the PDW stripes is the gapping of single-particle excitations 
in the antinodal region, as illustrated in Fig. IH in contrast, the nodal states remain 
ungapped. From the underlying band-structure, one sees that the largest contribution 
to the density of states with energies near Ep comes from the antinodal regions (where 
the dispersion is relatively flat); thus, the onset of local PDW order should have a major 
impact on properties sensitive to the total density of states. Conversely, properties that 
are largely determined by near nodal quasiparticle dynamics, which presumably includes 
the quasiparticle contribution to the in-plane conductivity, may be less strongly affected. 

Observed striking changes in various transport properties of several stripe order 
cuprates can be interpreted in this light as being suggestive of the appearance of local 
PDW order at the onset of charge-stripe order at Tco (which is generally somewhat 
higher than Tsdw)- In La2_a;Ba2;Cu04 and Nd- and Eu-doped La2_a;Sra,Cu04, it is 
observed that the in-plane thermopower drops dramatically below Tco [B 11201 11211 
I122j as does the Hall resistivity |123i 11241 I125j . Furthermore, the opening of a 
superconducting-like gap as the temperature drops below Tco results in an observed 
[126] suppression of the in-plane optical conductivity at frequencies below 40 meV. 
In contrast, the in-plane DC-resistivity changes relatively little [H 1127] upon cooling 
through Too- 

Putting aside the issue of the onset-temperature, the notion that stripe ordered 
cuprates exhibit local PDW order is also supported by ARPES studies. For example, 
measurements on stripe-ordered Lai.48Ndo.4Sro.i2Cu04 at T = 15 K (> 2Tc) reveal a 
gapless nodal arc of states covering roughly a third of the nominal Fermi surface, as 
well as a gap reaching 30 meV in the antinodal region [128J . Temperature- dependent 
ARPES measurements on La2_xBaxCu04 with a; = | indicate that, for temperatures 
above the spin-ordering transition, there is a gapless nodal arc of states, together with 
a substantial antinodal gap [15]. 

However, there are several aspects of this story which require further analysis. 
Firstly, there is the issue that different aspects of the crossovers we would like to identify 
with the onset of local PDW order appear to onset at different temperatures. This is not 
necessarily inconsistent, as a crossover (as opposed to a phase transition) can appear to 
occur at somewhat different temperatures depending on what quantity is measured and 
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how the data is analyzed. Nonetheless, the drop in the thermopower and Hall number 
appears to have a very sharp onset at Tco, while the superconducting like drop in the 
in-plane resistivity at Tsdw is also very sharp, at least in 1/8 doped La2-xBa^Cu04. 
[In this sense, it is reminiscent of the situation |129] in O doped La2Cu04, where the 
sharply defined spin ordering and superconducting ordering transitions occur at the 
same temperature (in zero field) with very small uncertainty.] 

A still more perplexing issue arises in correlating the onset of the signatures of 
2D superconductivity in Lai.875Bao.i25Cu04 with the thermal evolution of the ARPES 
(m Us] spectrum. Below Tsdw, there is clear evidence of the appearance of a d- 
wave-like gap in the nodal region, with the scale of this second gap being smaller 
than the pre-existing antinodal gap |15]. This behavior suggests that uniform ci-wave 
superconductivity develops and coexists with the PDW superconductor below Tsdw- 
However, this is somewhat problematic, as the proposed explanation of the dynamical 
interlayer decoupling and the bounded growth of superconducting correlations that 
occurs below Tsdw rests on the assumed (near) absence of a uniform component of the 
order parameter in each plane. Reconciling the uniform ci-wave component of the order 
parameter inferred spectroscopically from ARPES studies of La2-xBaa.Cu04 with the 
apparently almost complete absence of such a component inferred from bulk transport 
measurements on the same material is a challenge for future work. It may be significant, 
however, that ARPES studies of Lai.6-xNdo.4Sr3;Cu04 |128j and Lai.8_2;Euo.2Sra;Cu04 
|130j appear consistent with pure PDW order, {i.e., there is no d-wave gap in the nodal 
region), although the PDW appears to set in at around Tco, which can be substantially 
greater than Tsdw in these materials. 

9.2. Dynamical layer decoupling and quasi-two-dimensional behavior in the cuprates 

The cuprate superconductors are layered materials with varying degrees of quasi-two- 
dimensional behavior. Evidence for quasi-2D behavior (and for dimensional crossover) 
in the cuprates has existed for a long time and it is well documented. It is thus useful 
to compare and contrast this well known behavior with the unexpected layer decoupling 
effect observed in La2-a;Baa;Cu04. 

In a quasi-2D system, as a continuous thermodynamic superconducting phase 
transition is approached, the in-plane correlation length grows very rapidly. While 
at first the fiuctuations have a markedly 2D character, very close to the phase transition 
they rapidly cross over to their ultimate 3D behavior. Dimensional crossover is 
observed, for instance, in dynamical probes of some cuprates. High frequency (~100 
GHz) conductivity measurements in Bi2Sr2CaCu208+5 (the most quasi-2D material 
among the cuprates) by Corson et al |131j showed that (at those frequencies) the 
fiuctuation conductivity is 2D-like and exhibits Kosterlitz-Thouless behavior, as if the 
Cu02 planes were effectively decoupled. Similarly, quasi-2D behavior in the dynamic 
conductivity (with frequencies in the range 1-10 GHz) has been observed in underdoped 
La2_2^Sr2^Cu04 near Tc (but not in overdoped La2_a;Sra,Cu04) by Kitano et al |132j . By 
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probing the system at finite frequency, tliese experiments explore tlie correlations at 
a frequency dependent mesoscopic length scale, where sufficiently weak 3D couplings 
have negligible effect on the physics. By contrast, the resistive transition both in 
Bi2Sr2CaCu208+5 and in La2_a;Sra;Cu04, measured at zero frequency in macroscopic 
samples, is not of the 2D XY (Kosterlitz-Thouless) type, but rather reflects the three- 
dimensional nature of these materials. 

In contrast, the unusual layer decoupling effect observed in stripe- ordered 
La2-xBaa;Cu04 takes place in a temperature range where the Cu02 planes appear to 
become superconducting (well above the three-dimensional critical temperature). [H [28] 
The layer decoupling effect is observed in the resistive transition, and is thus not a 
dimensional crossover effect. As we noted above, in this regime La2_xBa^.Cu04 behaves 
as if for some reason the effective inter-layer Josephson coupling is either turned off 
(which is unphysical) or is somehow frustrated. 

Support for this idea is provided by recent Josephson resonance experiments in 
La2-xSr2;Cu04 by Schafgans et al[2], which essentially measure the c-axis superfluid 
stiffness, pc- In the absence of an external magnetic field, pc has the expected |133[ 
11341 1135] magnitude, i.e. pc is proportional to the normal state conductivity at Tc. 
However, for underdoped materials, pc becomes unmeasurably small in the presence of 
moderate magnetic fields {B < 8T). Magnetic fields are known to induce static spin- 
stripe order (as detected by neutron scattering experiments [H]) in precisely the same 
range of field strengths and hole concentration. These experiments thus suggest [2] that 
the "fluctuating stripe order" p!5] seen in La2-xSr2^.Cu04 at zero field may actually be 
of the PDW type and that dynamical layer decouphng occurs as static stripe order 
is stabilized in a magnetic field. § Indeed, in materials, including La2-xBa2;Cu04 and 
Lai.6-xNdo.4Sr2:Cu04, which exhibit stripe order in zero field, pc is found [T36] to be 
orders of magnitude smaller than its "expected" value on the basis of the normal state 
conductivity. 

9.3. Possible relevance to other cuprates 

Although there are still open issues, the PDW state (or its glassy version) seems to offer 
a rather compelling explanation for what is otherwise an extremely surprising set of 
phenomena observed in stripe ordered cuprates. Could these ideas also be relevant to a 
broader range of phenomena in the cuprates? The direct empirical information available 
[15] concerning the structure of any sort of static or fluctuating stripe order present in 
cuprates outside the 214 family is much less clear. || Consequently, any attempt to 

§ While it is tempting to reinterpret in hindsight the resuhs of Corson et al [131j as being indicative of 
"fluctuating PDW order" in Bi2Sr2CaCu208+5 we should note that the STM data on this material|5] 
show a glassy pattern of short range stripe order at high bias. However, as we explained elsewhere in 
this paper, a glassy version of the PDW state would not exhibit a sharp layer decoupling effect. 
II The results of recent neutron scattering studies of underdoped YBa2Cu306+a by Hinkov et al. [55] 
have confirmed |52j the existence of a nematic phase, derived from the weak melting of a stripe ordered 
state, onsetting below a temperature comparable to the pseudogap onset-temperature, T*. Still more 
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achieve a theoretical understanding based on the assumed existence of a PDW state is 
necessarily speculative. We therefore present the discussion of this final section in the 
spirit of provocative conjectures, which we believe are deserving of further investigation. 

ARPES studies of Bi2Sr2CaCu208+5 [82] and La2_xSr^Cu04 M have 

revealed "Fermi arcs" of gapless states between antinodal pseudogaps. There has been 
a great deal of controversy over the nature of the antinodal pseudogap |139[ 1140] . Two 
recent studies of Bi2Sr2CaCu208+5 have reported signatures of Bogoliubov quasiparticles 
in the antinodal gap region [63 lll41j . which was interpreted as being suggestive that the 
pseudogap is, at least in part, produced by superconducting fluctuations. On cooling 
through Tc, a rf-wave gap appears along the nodal arc |137l 11381 1142j . In near optimally 
doped samples, as T — > 0, this d-wave gap and the pseudo-gap merge to form a single 
gap with a simple [cos{kx) — cos{ky)] form. However, in underdoped samples, even as 
T — i> 0, the nodal gap appears to have a different energy scale than the antinodal gap 
{i. e., they do not merge to form a simple d-wave gap) [138U142] . Thus, in some ways it is 
clear that there are two distinct gaps - an antinodal pseudo-gap that might be associated 
with some sort of "competing" order, and a nodal gap, which is clearly superconducting 
in the sense that it onsets quite sharply at Tc. However, in other ways it seems that all 
the gaps have some unifying superconducting character. 

We propose that this puzzle may be resolved by postulating that there are two 
distinct gaps, both with superconducting character in the sense that one is associated 
with uniform the other with modulated superconducting order. Indeed, the measured 
quasiparticle spectral function in the pseudogap looks somewhat like that of the PDW 
state (see Fig. H]).^ Moreover, just such a combination of modulated and uniform 
superconducting orders has been previously proposed on phenomenological grounds 
to explain [62 [92] STM spectra ^ EM Mi [M [US] in Bi2Sr2CaCu208+5 and other 
cuprates [3 1146] . 

Seemingly more direct evidence of superconducting fluctuations in the normal 
state of La2_2:Sr^Cu04, Bi2Sr2CaCu208+5, and Bi2Sr2-j,Laj,Cu06 has been reported 
by Ong and coworkers [1471 1148] based on measurements of the Nernst effect and 
diamagnetism. Nernst measurements on YBa2Cu306+a; [149] and STM studies of 
Bi2Sr2CaCu208+5 [150] suggest that disorder may be important to the existence of 
the fluctuation effects over a substantial temperature range. It is intriguing that the 
onset temperatures of the enhanced Nernst response in La2_a;Sr2;Cu04 has a maximum 
at a; ~ 0.1 [147] . close to the optimum doping for stripe order. Moreover, Taillefer 

recently 156] . the same authors have demonstrated that modest magnetic fields stabilize static spin- 
stripe order where primarily fluctuating (nematic) order existed at zero field. 

% Technically, the electron-hole-mixed quasiparticles in the antinodal region of a PDW state are not 
perfectly symmetric with respect to Ep (see Fig. |4]), in contrast to Bogoliubov quasiparticles; however, 
to detect this distinction, one would need to measure a sample containing a single-domain PDW state. 
Measurements on a nematic PDW state in Bi2Sr2CaCu208+5 would average over stripe orientations; 
furthermore, the experimental "quasiparticle" peaks are quite broad in the pseudogap state [63 l 1141] . 
so that any fine details are hidden by damping. The overdamping also fills in the spectral weight at 
Ep, in contrast to the true gap that is found for T <Tc [63 11137] . 



Striped superconductors 



38 



and coworkers |151j have found close correlations between an enhanced Nernst signal 
and stripe order. Neither the observed sensitivity to disorder nor the association with 
stripe order, by themselves, necessarily negate the interpretation of these effects in 
terms of superconducting fluctuations; however, both would be unusual in the case of a 
simple, homogeneous d-wave superconductor. While we are far from having an explicit 
theory, it seems to us that these general trends are consistent with the existence of 
a disordered PDW state over at least a portion of the pseudogap phase. Specifically, 
Ong and coworkers [152] have reported the observation of a sublinear dependence of 
the magnetization on magnetic field (M ~ —5" with a < 1) in a relatively narrow but 
non-vanishing range of temperatures above Tc in crystals of Bi2Sr2CaCu208+5. This 
behavior, if it truly persists in the limit B 0, must signify the existence of a distinct 
phase of matter in this range of temperatures, which we very tentatively propose could 
be a superconducting glass formed from a disordered PDW.'^ 

One of the most intriguing recent discoveries in the cuprates involve several distinct 
observations of a rather subtle, and not fully understood, form of time reversal symmetry 
breaking in the pseudogap phase of YBa2Cu306+2; [571 158] HgBa2Cu04+5 [59]. As 
we have seen, various forms of subtle time-reversal symmetry breaking can occur when 
frustration is added into the PDW mix. It is our hope that, with further work, a relation 
can be established between these two rather vague statements. 
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